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Introduction 


While the role of exactly solvable and partially (quasi-exactly) solvable models in one¬ 
dimensional Quantum Mechanics is well known'^, the number of these models is rather re¬ 
stricted. For two-dimensional systems (exact- and partial-) solvability is a much harder 
task. Besides models with separation of variables, which are actually reduced to couples of 
one-dimensional models, only Calogero and the so called Calogero-like models[1] (with the 
number of particles N = 3)® are exactly solvable. Some partially solvable systems were con¬ 
structed in [3], [4], [5], [6] by the supersymmetrical method of S'f/S'F—separation variables. 
Just the methods of supersymmetry seem to be the most appropriate ones to attack this 
fundamental problem of constructing exactly and partially solvable quantum (and classical 
[7]) models in the case of few space dimensions. 

The hrst step in this direction was taken in [8], where a list of partial solutions of super- 
symmetrical intertwining relations for two-dimensional systems was found. All these models 
are integrable, i.e. there is a dynamical symmetry of fourth order in momenta, such that 
the corresponding generators are in involution with the Hamiltonians. After that two new 
supersymmetrical methods - of S'f/S'H-separation of variables and of shape invariance - were 
proposed [3], [5], [6], [9] which provided an opportunity to obtain part of the spectrum and 
wave functions for some of these integrable models. Thus a class of partially solvable inte¬ 
grable two-dimensional quantum systems was built. Naturally, it would be very important 
to obtain some exactly solvable models among these integrable ones. 

In the present paper a new idea, which may be useful in order to solve the task formulated 
above, will be presented. Namely, we will choose a particular value of parameter (a = —1/2) 
in the two-dimensional Morse potential with additional sinh“^ term [3] in such a way that one 
'^Here and below we call the system exactly solvable if all eigenvalues and eigenfunctions of its bound 
states are known analytically. The system is called partially (quasi-exactly) solvable when a part of them is 
known. 

^Incidentally, we would like to mention that like the Calogero models all considered models (including 
this generalized Morse model) cam alternatively be interpreted [2] either as describing three particles on a 
line or as a particle in a two-dimensional space. 
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partner Hamiltonian in the SUSY intertwining relations will allow for standard separation 
of variables. In such a case this Hamiltonian is exactly solvable, and its spectrum and 
eigenfunctions will be found analytically. The second partner Hamiltonian still does not 
allow for the separation of variables, but due to intertwining relations its spectrum and 
eigenfunctions can also be obtained analytically by acting with the supercharge. Since SUSY 
intertwining relates the spectra of partner Hamiltonians up to normalizable zero modes of 
supercharges, the requirement of exact solvability of both Hamiltonians amounts to the 
condition that all zero modes of supercharges are under control. This could be implemented 
for the generalized two-dimensional Morse model (with additional sinh“^ term), but as will 
be shown in Section 2, the direct analysis of singularities of wave functions for a = —1/2 
case is rather problematic in a general form due to the presence of confluent hypergeometric 
functions. Nevertheless the principal opportunities offered by the method are demonstrated 
for a particular energy level. Consideration in Section 3 of some specihc complex version 
of this model (which regularizes the repulsive singularity for the real values of xi,X 2 ) leads 
to a more interesting system: all bound states and normalizable wave functions can be 
derived explicitly. The properties of this last model are investigated. The most important 
feature is the reality of the energy spectrum (which does not depend on the parameter of 
complexihcation 5). In addition, in Section 4 the bound-state-biorthogonal-basis [10] will be 
constructed in a natural way, and the action of symmetry operator will be calculated. In 
particular, it will be shown that the obtained wave functions are the common eigenfunctions 
both of the Hamiltonian and of the symmetry operators. 

2. Generalized two-dimensional Morse potential with 

a = -1/2 

The main element of the two-dimensional second order Supersymmetrical Quantum 
Mechanics (SUSY QM) [8], [5] is represented in terms of the intertwining relations: 

HQ+ = Q+H; Q-H = HQ-- (1) 
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between a pair of scalar Hamiltonians H, H : 

H = -A + V{x); H = -A + ^(f); A = di = d/dxi] i = l,2. (2) 

These intertwining relations realize the isospectrality up to zero modes of of the super¬ 
partners H, H and the connection between their wave functions with the same values of 
energy: 

^n{,x) = Q~i)n{,x)] ^n{x) = Q^'^nix)-, n = 0,1,2,... (3) 

Though the solutions of intertwining relations (1) have to be searched in principle for the 
supercharges with the most general form of second order differential operators (see [8], 
[5] for details), we will restrict ourselves to the particular solution of Lorentz (hyperbolic) 
type: 

Q+ = {dl - di) + CA + B = 4d+d_ + C+d_ + C_d++ B- (4) 

Q- = {di - dl) - Cidi + B = 4d+d_ - C+d_ - C_d+ + B- i = l,2. (5) 

Then the potentials V{x),V{x) and the function B{x) can be expressed in terms of four 

functions - Fi(2a;i), ^ 2 ( 2 x 2 ) and C±{x±) : 

V = i(c; + Gi) + l(c|+ci) + i(f;,(i+-i_)-f’i(2^+ + 2^-)). 

r = -i(C;+C^) + i(Ci + Gi) + l(F2(i+-i_)-Fi(i+ + i_)). (6) 

B = t^C+C_ + Fi(x++x_) + F2 (x+- x_)). (7) 

These functions must satisfy the following equation: 

d.{C.F) = -dAC+F), 

where x±=xiA X 2 d± = d/dx± and C± depend only on x±, respectively: 

C'+ = Ai - ^2 = C+{xA C_ = Ci + C2 = C_{x-), 

and 

F = Fi{x+ + X-) + F2 {x+ - X-). 
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Among different particular solutions of this class [8], [5] we will choose here the well 
studied two-dimensional generalization of one dimensional Morse potential. The model is 
not amenable to standard separation of variables, and it is dehned by: 



= 

(y.x 

4aa] C_ = 4aQ; ■ coth 


(8) 


= 

1 / ^ 

-Fi(2a;i) = — A( exp(—2Q;a;i) — 2exp{—axi) 
4 \ y 

h 

(9) 

f2ix2) 

= 

1 / ^ 

-^ 2 ( 2 x 2 ) = +A[ exp{-2ax2) - 2exp{-ax2) 
4 \ y 

) 

(10) 

V{x) 

= 

a‘^a{2a — 1) sinh“^^^^^^ -|- 4a^Q;^ -|- 




+ 

A exp(—2Q;a;i) — 2exp(—(axi) -|-exp(—2Q;a;2) 

— 2 exp(—Q;a; 2 ) 

(11) 

V{x) 

= 

a^a{2a + 1) sinh“^ 




+ 

A exp(—2Q;a;i) — 2exp(—axi) + exp(—2ax2) 

— 2 exp(—Q;a; 2 ) 

(12) 


where parameters a, A > 0, a > 0 are arbitrary real numbers. 

Just on the basis of this model two new methods - S'l/S'K-separation of variables and 
two-dimensional shape invariance - were elaborated for the first time [3], [5]. As a result, 
partial solvability of the model was discovered [3], [5], [9] and a variety of wave functions 
was found analitycally for parameters restricted to the range: 

. 1 1 , VA 1 

a E (—oo, - -^);- n -> —2a >0; n = 0,1, 2,..., (13) 

4 4^2 a 2 ^ ^ 

which provides the condition of normalizability of zero modes of Q~^ and the absence of the 
’’fall to the centre” (see details in [3],[5]). 

It is worth to recall now the main idea of the method of S't/AF-separation of variables. 
By similarity transformation the supercharge (4) for Lorentz metrics can be transformed 
to the operators without linear derivatives: 


= exp(-x(T))Q^exp(+x(T)) = (9^ - (92 + ^(Fi(2a;i)+^ 2 ( 2 x 2 )); (14) 

X(x) = -^(^J C+{x+)dx++ J C_{x_)dx_y 
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which allow for the separation of variables (this is why we call this method as SUSY- 
separation). If —Fi and +F 2 belong to a class of exactly solvable one-dimensional potentials, 
the normalizable zero modes of Q~^ can be found analytically. Then due to (1) a set of 

wave functions of H can be constructed as linear combinations of these zero modes, leading 
therefore to partial solvability of the model. The shape invariance method [3] (together 
with the second shape invariance [9] of the model) enlarges the part of the spectrum known 
analytically. 

Now a new method, involving separation of variables in one partner Hamiltonian (which 
does not hold for the other partner), will be formulated. Let us choose the parameters of the 
same model in such a way that the Hamiltonian H does allow for the standard procedure 
of separation of variables. Then the intertwining relations and knowledge of zero modes of 
supercharges provides full information about the partner Hamiltonian H, which does not 
allow conventional separation of variables. 

Luckily, a suitable choice of the value of the parameter a = —1/2 in (12) makes H 
amenable to separation of variables: 

H{x) = hi{xi)+ h 2 {x 2 ) + hi{xi) =-dl - fi{xi)] h 2 {x 2 ) =-dl + f 2 {x 2 ), (15) 

and its wave functions with energies 

En,m = Y €m O? ■, (16) 

where the last term originates from the free term in (12), can be written as: 

= Cir]n{xi)r]m{x2) + C2r]m{xi)r]n{x2)] Cl, C 2 = CoUSt (17) 

where and rik{x) solve the exactly solvable one-dimensional problem with the standard 
Morse potential in terms of confluent hypergeo metric functions (see [11]): 

^—(9^ -|- H^exp(— 200 ;) — 2exp(—Q;a;)^^?7fc(a^) = ^kVk{x) (18) 

c 2\/A 

% = exp(-|)(0"''$(-/c,2sfc + 1;0; ^ =-exp(-aa;); (19) 

2 a 

ek = -A[l-^{k + l/2)]^; s^ = ^-k-l/2. ( 20 ) 

y A o> 
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Thus the two-dimensional Schrodinger problem with Hamiltonian H is obviously exactly 
solvable with 2-fold degeneracy for the levels with n ^ m. 

Now the next step of our approach is to use the intertwining relations (1) in order to 
obtain, by means of (3), the wave functions of the Hamiltonian H, which does not allow for 
separation of variables. It is clear from the explicit form of Q~^ that due to the singularity of 
Q~^ at a;_ —0 the normalizability of functions T should be investigated in detail. Though 
this problem is not analyzable fully due to the presence of hypergeometric functions in (17), 
particular normalizable eigenfunctions of H can be studied. Indeed, taking antisymmetric 
(ci = —C 2 ) wave functions of the form (17), one can check^ that for |n — m| = 1 they 

are represented as linear combinations of zero modes VLk{x) of Q^\ 




2(<Sm+l + l)(2Sm+l + 1) n / 

2s^^^ + m + 2 


and therefore have no partner bound states at H. The symmetric functions (ci = C 2 ) from 
(17) with n = m±l give which are also absent among bound states of H due 

to their nonnormalizable behaviour at X- ^ 0. Therefore the spectrum of H definitely does 
not include Em±i,m- 

The wave function leads to the (symmetrical) which can be rewritten as: 

'i'l.,, = = “'(2*2 + 3)K. - exp(-h^)K,^j)». 


This example demonstrates explicitly that a suitable choice of ci,C 2 ,n,m at (17) can com¬ 
pensate the singularity in Q~^ and provide a normalizable wave function 

In principle, besides the eigenfunctions T of 77 discussed above, some additional nor¬ 
malizable eigenfunctions could exist. If so, due to intertwining relations the action of the 
operator Q~ onto these functions should either give zeros or unnormalizable functions T. 
The first option is excluded by the analysis of zero modes of Q~{a = —1/2) = Q~^{a = -1-1/2), 
taking into account that a = -1-1/2 is outside the range (13). The second option also does 

tit is necessary to use elementary relations between confluent hypergeometric functions (see [11], vol.l, 
Subsection 6.3.) and some properties derived in Subsection 4.4. of [3]. 
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not materialize since normalizable eigenfunction of H with potential (11) for a = —1/2 gives 
tl' ~ at a;_ —0. Therefore action of Q~ cannot affect its normalizability, and we are 
back to the case of previous paragraph. 

3. Complexified model 

Since the bound state spectrum of the model i/ for a = —1/2 turned out to be difficult to 
analyze due to singularities of V{x) and at a;_ = 0, a natural idea is to remove somehow 
these singularities. In one-dimensional Quantum Mechanics the recipe is well known [12], 
[13], [14] - to shift the space coordinate into the complex plane. It means that from this 
moment we deal with complex potentials^, in general. 

In our two-dimensional situation^ we have to shift X- = xi — X 2 , therefore one has to 
violate the exchange symmetry of the system under xi ^ X 2 - The easiest way is to replace 

X ^ X + id; 5 = {5, 0) 

(with S small enough, such that aS G (0,7r/2)) removing the singularities from the real 
{xi,X 2 ) plane. In terms of ^ of Eq.(19) a phase factor appears: ^ As usual, such 

imaginary shift of x preserves the reality of the spectrum of the Schrodinger operator. 

Under this shift not too many changes affect formulas of Section 2. The operator Q~{x-\- 
iS) in (5) preserves its form, but from now on it is not hermitian conjugate of Q~^{{x -1- id)), 
namely Q~ = = (Q"*")*- Functions (8) - (10) become complex, and (also complex) 

potentials (11), (12) are related now by 

V{x + id) = V*{x — id) = exp {—2idd)V*[x + id) exp {+2idd). 

This equality expresses the property of pseudo-Hermiticity for non-Hermitian Hamiltonians, 
which guarantees in general [10] that the spectrum consists of real eigenvalues and complex 
®An extensive literature concerning one-dimensional non-Hermitian Hamiltonians followed the seminal 
papers [15] of C.Bender and S.Boettcher (see for example, [16]). 

^Non-Hermitian models in two-dimensional Quantum Mechanics were studied in [4]. 



conjugated pairs. In particular, in our case the whole bound states spectrum of H is known 
to be real. 

The eigenfunctions r]k{x) of one-dimensional Morse equation (17) are expressed in terms of 
confluent hypergeometric functions as in (19), but with ^ and are still normalizable. 

One can check that no additional normalizable solution of (17) appears in the complex x- 
plane. Indeed, after the substitution rjk = exp(—,^/2),^®y(,^), the Eq.(18) in the variable 
^ is reduced to the confluent hypergeometrical equation [11]. There are different ways to 
represent the general solution of this equation. A convenient one is: 

= Cl|/5(0 +C2|/7(0, 

where ^ includes a phase factor, and dehnitions and terminology for linearly independent 
solutions 1 / 5 , 1/7 are given in [11] (see Vol. 1, Subsection 6.7.): 

//5 = $(a,&;0; 1/7 = exp(0T(6-a,6;-0- (21) 

Just the exponential in the second term in (21) allows to prove that the only kind of normal¬ 
izable solutions F, even for our complex corresponds to C 2 = 0. Therefore, the conditions 
of normalizability lead to Eqs.(19), (20), and the whole bound state spectrum (16) of Hamil¬ 
tonian H{x -\- id) is known and it remains real after this complexihcation. 

Using again the intertwining relations (1), one can obtain the eigenfunctions 

+ id) = Q^r]n{xi -f id)7]m{x-2). 

Due to absence of singularity of at a;_ —0, these wave functions are normalizable. 

One has to remember that the partner Hamiltonians H and H are isospectral only up 
to zero modes of supercharges. Normalizable zero modes of were constructed and in¬ 
vestigated in detail in [3]. It was proved that a variety of linear combinations of these zero 
modes can be built, which are eigenfunctions of the Hamiltonian H. In particular, this can be 
done for the case a = —1/2 (when H allows separation of variables), and the corresponding 
eigenvalues are [3]: 

\fA 

Ek = 2efc -|- 2a^Sk = —2a^Sk{sk — 1) = e^+i -f- J- Sk =- k — 1/2. (22) 

a 
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For arbitrary n,m eigenvalues (16) of H are two-fold degenerate: there are symmetrical 
and antisymmetrical components in (17). But for the particular case of n = m ± 1 the 
antisymmetrical combination annihilated by . In a contrast to the Hermitian 

model of Section 2, for the symmetrical component its partner state = 

the complex model has no singularity at a;_ —0, and therefore these energy 
levels Em±i^m exist in the spectrum of H, but they are not degenerate. 

The possible zero modes of Q~{a) = Q~^{—a) can be also investigated by the same method 
as in [3], the corresponding eigenvalues are = —2a^Sk + 26^. Because these eigenvalues 
coincide with i7fc,A:-i from the general expression (16), no new eigenstates can appear for H 
due to zero modes of Q~{a). 

Thus the complex model with the Hamiltonian H{x + i6) for a = —1/2 is exactly 
solvable, its spectrum is real: 

En,m = Cn + (23) 

For n = mil these levels are not degenerate, but for all other n, m there is 2-fold degeneracy. 
The wave functions of H are: 

(24) 

where E„,m sue given by (17) with x ^ x + iS; S = {6, 0). 

4. Integrability and biorthogonality 

It is known [8], [5] that ”by construction” all Hamiltonians 77,77, which are intertwined 
according to (1) by second order operators Q^, are in involution with operators of fourth 
order in derivatives: 

R = Q+Q--, 77 = g-g+; [77,77] = 0; [H,R] = 0. (25) 

These symmetry operators* are not reducible to functions of Hamiltonians, and therefore all 
these systems, including the one presented in this paper, are integrable. 

*Apart from the case of Laplacian metrics in supercharges, which is not considered here (see details in 

[ 8 ]). 
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In the particular case a = —1/2 with separation of variables in H, considered in this paper, 
the expressions for R and R can be transformed by using the similarity transformation (14) 
and the specihc form of one-dimensional Hamiltonians hi{xi), h 2 {x 2 ) in (15): 

= exp(±x)g^exp(=Fx) = exp(±x)(/i 2 - hi) exp(=Fx). 

Then the symmetry operator R for the Hamiltonian with separation reads: 

R = (^2 — hi)^ 2a‘^[hi ^ 2 ) + 

and its wave functions (17) are simultaneously eigenfunctions of R with eigenvalues: 

f'n,m = — ^n)^ + 2Q;^(em + ^n) + ■ 

The wave functions (24) of the Hamiltonian H without separation of variables are also 
common eigenfunctions both of H and R : 


Thus the property, noticed in [3] for a few known wave functions, is fulhlled now for all wave 
functions in the case of the exactly solvable model, that we have constructed in this paper. 
Though almost all (for n ^ m) wave functions '^E„m 2-fold degenerate, the symmetry 
operator R does not mix these degenerate functions. 

The factorized wave functions of 77(1; -f iS) and their complex conjugate 

functions '^Enm fh® so called biorthogonal basis for the non-Hermitian Hamiltonian. 
The corresponding biorthogonality relations 


< I >= J + i5) ■ + id) = 

= J dxir]n{xi + iS) ■ r]n'{xi + iS) J dX 2 Vn{x 2 ) ■ Vn'{x 2 ) = Snn'dn 


can either be checked straightforwardly or by comparing this integral along the line xi + id 
with the analogous integral along the real Xi for the case of <5 = 0 and real-valued wave 
functions ? 7 „. Absence of singularities in the narrow strip between these lines means that 
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relations (26) follow from the standard orthogonality of wave fnnctions for the Hermitian 
Hamiltonians. On the contrary, if one wonld consider the scalar prodnct dehned with the 
complex conjngation of the hrst mnltiplier in the integral in (26), the connection with the 
Hermitian case would be not so evident. 

The bound-state-biorthogonal-basis for the non-Hermitian operator H{x + iS) consists of 
and = (Q+)*d'^^^. Due to the equality Q" = ((Q+)^)* the scalar 

products can be written as: 


< T) 




>=< (g+)*T) 




>=< n. 


g-g+T 




Since is an eigenfunction of the symmetry operator R = g^g"*" with an eigenvalue 

rn,m, biorthogonality for H{x+iS) follows directly from (26), thus leading to a diagonalization 
of H in the bound state subspace. 
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